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Abstract 

Chimera is a relatively new emerging phenomenon where coexistence of synchronous and asyn¬ 
chronous state is observed in symmetrically coupled dynamical units. We report observation of the 
chimera state in multiplex networks where individual layer is represented by 1-d lattice with non-local 
interactions. While, multiplexing does not change the type of the chimera state and retains the multi¬ 
chimera state displayed by the isolated networks, it changes the regions of the incoherence. We inves¬ 
tigate emergence of coherent-incoherent bifurcation upon varying the control parameters, namely, the 
coupling strength and the network size. Additionally, we investigate the effect of initial condition on 
the dynamics of the chimera state. Using a measure based on the differences between the neighboring 
nodes which distinguishes smooth and non-smooth spatial profile, we find the critical coupling strength 
for the transition to the chimera state. Observing chimera in a multiplex network with one to one inter 
layer coupling is important to gain insight to many real world complex systems which inherently posses 
multilayer architecture. 


1 Introduction 


In past few decades, network science has discovered a plethora of novel phenomena while trying to mimic 
real world systems in a better manner. One of such discovery is an observation of the chimera state. It 
was first reported by Kuramato et. al. in 2002 while investigating non locally coupled identical oscillators 
in a ring network [Kuramoto & Battogtokh (2002) ]. Later, it was analyzed and christened by Abrams and 
Strogatz in 2006 as chimera state [Abrams & Strogatz (2006)] . Like, its counterpart in Greek mythology, 
a chimera state has come to be referred as a mathematical hybrid state in which coherent and incoherent 
dynamics coexist in non-locally coupled identical oscillators in a structurally symmetric network. 

Chimera has been extensively investigated both theoretically [Sethia et aL (2008)||Laing (2009) ,Omelchenko et_ al. 
and experimentally [Hagerstrom gf aL (2012)][Larger et al. (2013)] . It has been observed in plenty of net¬ 
works including phase oscillators [ Abrams et_ aL (2 004)), Abrams & Strogatz (2006)[Maistrenko et aL (2014)1 
Omel’chenko et al (2012)[Sethia et al. (2008)|Laing (2009)| , chemical [ Tinsley et aL (2012)||Nkomo et aL (2013)B , 
mechanical oscillators [Martens etaL (20 0)1 , neuron models [Hizanidis et al (2014)[ , planar oscillators 
[Laing (2010)] ], boolean networks [Rosin et aL (2014)||Rosin (2015)[ , ID superconducting meta material 
[Lazarides etaL (2015)[ , etc. Chimera was originally reported for non-locally coupled oscillators, recently 
it has been reported in feed back delayed networks [Omerchenko et aL (2008)[|Sheeba et al. (2010)] ], glob¬ 
ally coupled networks [Schmidt et aL (2014)[Yeldesbay et aL (2014)[ , time varying networks [Buscarino etaL (2015)] 
and networks with purely local coupling [Laing (2015)| . Moreover, different types of chimera has been re¬ 
ported including multi cluster chimera [Xie et al. (2014),Omelchenko et al. ( 2013)[ , virtual [Lar ger et al. (2 013)] ], 
breathing [Abrams et aL (2008)[ and two dimensional chimera [Omel’chenko et al (2012)[Panaggio & Abrams (2015 
A recent work suggests emergence of Chimera, dependent on nonhyperbolicity of dynamical systems for 
both the time-discrete and time continuous cases [Semenova et al. (2015)]. 
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Figure 1: Schematic diagram for multiplex network consisting of two layers. Each layer is represented by 
ID lattice with non-local interaction. Each node (circle) has the same coupling architecture. 


The chimera state has also been reported for various real world networks models such as Rosenzweig- 
MacArthur oscillators for ecological networks [jDutta & Bane rjee ( 2015)]. Chimera has also been charac¬ 
terized by the state of the dynamical evolution of the network. Type -1 chimera is characteristic of the hyper 
chaotic behavior with many positive Lyapunov exponents HWolfrum et ai (201 l)fl . This type of chimera 
has primarily been observed for time-continuous systems like complex Ginzburg-Landau oscillators or Ku- 
ramoto oscillators. In Type-II chimera, only spatial chaos has been observed with a rather simple temporal 
behavior (mostly periodic). Though, this type of chimera has been reported mainly for the time-discrete 
systems (maps) flOmelchenko et aL (201 l)[[Omelchenko et al. (2012)[|Hagerstrom et al. (2012)] , it has re¬ 
cently been observed for the time-continuous Stuart-Landau oscillators as well HZakharova etai (2014)| . 

Further, modeling real world complex systems under the multiplex network framework is one of the 
recent advancements in the network theory Poccaletti et al. (2014)[|Lee et al. (20l2)l|Wang et aL (201 5)[ 
Kivela etaL (201 4)1 . We consider a multiplex network consisting of the same nodes across the layers 


(Fig. H|) and investigate the occurrence of chimera state in the multiplex ring networks. In this structure, 
each node has exactly the same connection architecture. We observe Type-II chimera with spatial chaos 
and periodic temporal behavior. Though, the chimera state, upon multiplexing, remains of the same type 
as observed for the isolated network, the multiplexing changes the region of incoherence. Dependence of 
the chimera state on initial conditions is observed for the multiplex networks as already observed for the 
isolated networks. Additionally we present a measure in terms of distance variable to distinguish between 
the coherent and the chimera state and to identify the transition point for the coherent-incoherent bifurcation. 
We also investigate role of the size of the network in determining the critical coupling strength for the 
symmetry breaking and thus emergence of the chimera state. 


2 Model 


We consider a multiplex S 1 ring network with N nodes in each layer, where S 1 represents one dimensional 
symmetric cyclic group with elements being invariant to the permutation operation flJacobson (2009)| |. Con¬ 
sidering zt(i),i = 1,..., rriN as a real dynamical variable at time t for the i th node, the dynamics of the 
network can be described as, 


mN 

z t+ 1(*) = f(zt(i)) + (2WV + l) [/(**(•?')) - /(**(*))] (!) 

where £ represents the coupling strength, m represent number of layers, r represents the coupling radius 
defined by r = P/N, with P signifying the number of neighbours in each direction in a layer. The elements 
of the adjacency matrix A of a network is defined as A t] = 1 or 0 depending upon whether i th and j th 
nodes are connected or not. The diagonal entries A t] = 0 represents no self connection in the network. The 
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adjacency matrix A for the multiplex network can be written as, 


A 1 I . 
I A 2 I 
. I A 3 


I\ 
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where A l ,A 2 ,...,A m represent the adjacency matrix of the first, second,... ,m th layer, respectively and / is 
an unit NxN matrix. Note that the number of nodes in each layer of the multiplex network is same. A 
mismatch in the network size of the layers will yield nodes in different layers having different interaction 
patterns and hence we can not define the chimera state. We use logistic map f(z ) = fiz(z — 1) with the 
bifurcation parameter // = 3.8 at which individual logistic map exhibits chaotic behavior. We consider 
coupling radius r = 0.32 for each layer indicating degree of each node being 64. A state of the network 
is defined as spatially coherent if for any node i E S 1 , the spatial distance between the neighboring nodes 
approaches to zero for t —» oo 


lim | z t (i + 1) — Zf(i) |-» 0, V i G S 1 (2) 

t—> OO 

Geometrically, this signifies a smooth profile of the spatial curve. Smoothness of the curve, signifying the 
correlated spatial values of the neighboring nodes, is defined with the absence of any discontinuity in the 
spatial curve. Whereas, temporal coherence (synchronization) is defined as, lim | z t (j) — z t (i) | —» 0 for 

t—¥ oo 

V i, j G S' 1 . Therefore, temporal coherence can be written as z t (1) = z t ( 2) = . = z t which leads to a 

straight line for the spatial curve associated with the temporal coherence. Appearance of discontinuity in 
the smooth spatial curve implies coexistence of the coherence and incoherence regions. To demonstrate the 
absence of smoothness, we define a measure based on the spatial distance as follows, 

d(i) = ( z(i + 1) - z(i)) - ( z(i ) - z(i - 1)) (3) 

which captures the difference of the spatial distances between the neighboring nodes. For smooth spatial 
profile d(i) = 0, signifying a symmetric distribution of the distances between neighboring nodes, while 
discontinuity in the spatial profile, signifying the transition point, is indicated as kinks in the distribution. 
We use this measure to find the critical coupling strength for the symmetry breaking and thus resulting the 
chimera state. 


3 Coherent-Incoherent bifurcation 


We evolve Eq.Q] starting with a set of special initial conditions and after an initial transient, study the spatio 
temporal patterns of the multiplex network. Note that uniform or the Gaussian distributed random initial 
conditions lead to either a completely coherent, spanning all the nodes, or a completely incoherent state 
depending upon the coupling strength. For 0.28 < e < 0.32 leads to the incoherent evolution of all the 
nodes. Motivated from [ |Kuramoto & Battogtokh (2002)HAbrams & Strogatz (2006)| |, we use a hump back 
function to generate initial conditions as follows. We choose an uniform random number z 0 (i) for initial 
state for i th node within some interval which varies like a Gaussian function as: 


z 0 (i) = exp[ 


0 ^_ 2 _ 

2cr 2 


N \2 


-] 


(4) 
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Figure 2: Snapshots and spatio-temporal plots for multiplex S' 1 ring network, (a) and (f) are for e = 0.1, 
(b) and (g) e = 0.28, (c) and (h) e = 0.30, (d) and (i) e = 0.40, (e) and (j) e = 0.44. Number of nodes in 
each layer remains N = 100 and coupling radius r = 0.32. Initial transient is taken as 5000. results are 
presented for time range 5000 to 5015. 


The variance a is chosen depending upon the size of the network such that the random variable lies between 
0 and 1. The same initial condition is used for both the layers in the multiplex network. A very narrow 
width of the function leads to almost very close value of the initial conditions for a fraction of nodes leading 
to the coherent state. 

In the absence of any coupling between the nodes (e = 0) or for weak couplings, all the nodes evolve 
independently and no spatial coherence is observed. For instance, as demonstrated in Fig.|2ta), for e = 0.1 
, the evolution of the nodes in the multiplex network yields an incoherent state with no correlations in the 
neighboring nodes. As the coupling strength is increased, a partially coherent state emerges at e = 0.28 
with correlated spatial values of the neighboring nodes in the end and in the middle regions of each layer, 
however, the spatial range of the incoherent region is more than the coherent region (Fig. ®b)). This 
coexistence of the coherent and incoherent dynamics corresponds to the chimera state in the multiplex 
network. 

The dynamical behaviour of two layers of a multiplex network is a replica of each other manifesting ex¬ 
actly the same spatio-temporal patterns (Fig.®. Exactly same behaviour is observed for multiplex networks 
having more than two layers (Fig®). 

At the same coupling value, the spatio-temporal dynamics (Fig. [2] (g)) reflects non-regular skeletal type 
pattern in the incoherent regions. This irregularity of the pattern suggests that, in the multiplex network 
framework, a node may get attracted to either of the upper or lower region depending on its initial value as 
reported for the isolated network [ Omelchenko et_ al. (2011)]. 

As we increase the coupling strength further, range of the incoherent region decreases as depicted by 
Fig.®c) for e = 0.3. At e = 0.4, we observe a sharp discontinuity in the otherwise smooth profile of z(j) 
and the incoherency appears at two distinct points in each layer. This is a bifurcation point for the coherent- 
incoherent transition. Above this coupling value, all the nodes in the multiplex network acquire the complete 
coherent state as indicated by the appearance of a smooth geometric profile at e = 0.44 (Fig. ®e)). Fig® 
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Figure 3: Distance measure for multiplex S 1 ring network (a),(c) for e = 0.44 and (b),(d) for e = 0.4. 
Network size is N = 100 



Figure 4: Snapshots for (a) Three layer and (b) Four layer multiplex network. Parameters are e = 0.28 and 
r = 0.32. Number of nodes N = 100 for each layer. 


depicts spatial regions of incoherent nodes and thus indicates a non-zero spatial entropy with the periodic 
temporal dynamics, representing a Type II chimera state. Further, the regions of incoherence in the spatial 
profile continues to exist for narrower intervals with an increase in the coupling strength (Fig.©. 

Furthermore, in the Chimera state, the time evolution of all the nodes in the network depict periodic 
behaviour with the periodicity two depicting temporal regularity. The coupled dynamics displays the spatial 
chaos which is defined by the non-zero spatial entropy given by h = d log e (2), where d represents fraction 
of the incoherent nodes in network [ [Omelchenko et aL (201 l)[[Coullet et al. (1987)] . We show that the 
distance measure (Eq© is able to easily distinguish between coherent and chimera state. The discontinuous 
spatial profile (Fig. (Mb)) at e = 0.4 gives rise to the kinks (Fig. ©d)) in the distance measure distribution 
signifying transition to the chimera state. We calculate the spatial entropy as a function of the coupling 
strength (e) in order to demonstrate the transition between the chimera to the coherent state. A transition 
from the chimera to coherent state is indicated by the discontinuous change in spatial entropy of the network 
(Fig©. 


4 Multiplex network versus isolated network 

In order to see the impact of multiplexcity on the dynamical behavior of nodes, we compare the dynamical 
state of an isolated ID lattice with that of the multiplexed with another ID lattice. We find that while the 
chimera state is retained after multiplexing, the dynamical evolution of the network differs. The multiplex- 
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Figure 5: Spatial entropy as a function of e. Figure are plotted for N = 100 and r = 0.32 
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Figure 6: Snapshot and Spatio temporal Pattern for (a,c) Single and (b,d) first S' 1 ring of multiplexed network 
are shown. Parameters are e = 0.30 and r = 0.32. Number of nodes N = 100 for each layer. Initial transient 
is taken as 5000. results are presented for time range 5000 to 5015. 


ing may enhance or suppress the incoherency. For instance, at e = 0.30, the isolated network displays 
the chimera state with incoherence in the middle regions of the spatio temporal pattern (Fig. 0 (a) and (c)). 
After multiplexing, the region of incoherence shrinks to a point discontinuity in the middle and intermediate 
regions (Fig.©. Thus, multiplexing here retains the chimera state as well as the type of the chimera state, 
but leads to a change in the region of incoherence as well as in the dynamical evolution. 


5 Sensitivity to initial conditions 

Furthermore, similar to the isolated networks, in the multiplex networks as well, the chimera exhibits de¬ 
pendency on the initial conditions. As already discussed in the section III, the spatial profile displayed by 
the chimera state can be very different for different profile of initial conditions for the same set of control 
parameters, namely /i and r. But interestingly, even for the initial condition given by the same profile as 
Eq. 0] with a constant value a for a given network size N, different realizations of the initial conditions 
can lead to a different incoherent regions. For example, at s = 0.30, for three different realizations of 
the initial conditions, all given by Eq. ffi three different spatio temporal patterns are observed (Fig. [7]). 
Though, the multi-chimera state is evident for all the realizations, the region of incoherence differs without 
any consistent behavior. 

Moreover, we investigated the impact of network size on emergence of chimera state as well as on the 
critical coupling strength below which one observes incoherent-coherent regime. Fig. [8] presents chimera 
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Figure 7: Spatio temporal pattern for multiplex S 1 rings for different realizations of the initial conditions. 
Parameters are e = 0.32 and r = 0.32. Number of nodes N = 100 for each layer. Initial transient is taken 
as 5000. results are presented for time range 5000 to 5015. 
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Figure 8: Spatio temporal pattern for multiplex S 1 rings with coupling strength (a) e = 0.21 (b) s = 0.23 
(c) e = 0.28 and (d) £ = 0.27 . Coupling range r = 0.32 and network size for (a),(c) N = 500 and (b),(d) 
N = 1000. Initial transient is taken as 5000. results are presented for time range 5000 to 5015. 


state for a multiplex ring network with two different network size. For both the network sizes, the cou¬ 
pled evolution exhibits the coexistence of coherence-incoherence dynamics. However, the critical coupling 
strength for the coherent-incoherent bifurcation increases to £ = 0.41 (Fig® as compared to £ = 0.4 for 
network size N = 200 as indicated by Fig.[2l 


6 Conclusion 

To summarize, we report an emergence of the chimera in the multiplex networks with the layers being rep¬ 
resented by 1-d lattice architecture having non-local couplings. We find that an emergence of the chimera 
is identical in the mirror layers arising due to the underlying symmetry of the network. Furthermore, while 
the temporal behavior of the network remains periodic even after multiplexing, the range of the coupling 
strength for which chimera is observed changes. The chimera in the multiplex network is found to be 
sensitive to the changes in the initial conditions, which is revealed through the changes in the incoherent 
region of the dynamical evolution for different sets of initial conditions. We also show that the critical 
coupling strength increases with the size of the network. The results presented here may provide a bet¬ 
ter understanding to the peculiar nature of the chimera state observed in many natural systems like uni- 
hemispheric sleep, ventricular fibrillation, brain networks which incorporates multilayer network architec- 
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Figure 9: Distance measure for multiplex S 1 ring network. (a),(c) for s = 0.44 and (b),(d) for e = 0.41. 
Network size N = 500 and coupling range r = 0.32 in each layer. 


ture [Panaggio & Abrams_(2015j|. 
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